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ABSTRACT 

We demonstrate that gamma-ray burst afterglow spectra and light curves 
can be calculated for arbitrary explosion and radiation parameters by scaling 
the peak flux and the critical frequencies connecting different spectral regimes. 
Only one baseline calculation needs to be done for each jet opening angle and 
observer angle. These calculations are done numerically using high-resolution 
relativistic hydrodynamical afterglow blast wave simulations which include the 
two-dimensional dynamical features of expanding and decelerating afterglow blast 
waves. Any light curve can then be generated by applying scaling relations to the 
baseline calculations. As a result, it is now possible to fully fit for the shape of the 
jet break, e.g. at early time X-ray and optical frequencies. In addition, late-time 
radio calorimetry can be improved since the general shape of the transition into 
the Sedov- Taylor regime is now known for arbitrary explosion parameters so the 
exact moment when the Sedov- Taylor asymptote is reached in the light curve is 
no longer relevant. When calculating the baselines, we find that the synchrotron 
critical frequency Um and the cooling break frequency u^. are strongly affected 
by the jet break. The Um temporal slope quickly drops to the steep late time 
Sedov- Taylor slope, while the cooling break z/c first steepens then rises to meet 
the level of its shallow late time asymptote. 



Introduction 



Gamma-ray bursts (GRBs) are currently thou ght to result from the collapse of a massive 



Eichler et al 



star ( Woosleylll993l : iMacFadven fc Wooslev 1999) or a black hole-neutron star or neutron 



1989 



Paczynskil 1 199 if ) . During these processes. 



star-neutron star merger (e.g 
a coUimated relativistic blast wave is launched into the circumburst medium. The emission 
from the blast wave is commonly referred to as the afterglow of the burst and can be ob- 
served throughout the broadband spectrum as the blast wave decelerates and radiates at 
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progressively longer wavelengths ( iMeszaros &: Reed 119971 ) . Ever since the discovery of the 
first afterglows, these have been modeled succesfully by combining a model for the blast wave 
dynamics with a sychrotron radiation model, where s hock-accelerated particles radiate by 



inter a cting with a shock-generated magnetic field ( e.g. iWijers et al.lll997l : IWijers fc Galama 
19991 : iFrail et al.l l2000t iPanaitescu fc Kuman |2002| ) . Analytically tractable solu tions for the 
dynamics are the self-simil ar Blandford-McKee (BMjBlandford fc McKedll976l ) and Sedov- 



von Neumann- Taylor (ST, ISedovlll959l : lTayloiill950l : IVon Neumann 



19611 ) solutions describ- 



ing, respectively, the ultra-relativistic and non-relativistic phase of the blast wave evolution. 
At early time lateral spreading of the collimated jet has not yet set in and the outfiow is 
purely radial, while at late times the jet will have become truly spherical, allowing the ap- 
plication of spherically symmetric solutions in both cases. As of yet, no analytical solution 
exists that fully captures the intermediate stage of the blast wave evol ution, where the blast 
wave becomes transrelativistic, inhom ogenous along the shock front (jZhang &: MacFadyen 
20091 : 1 Van Eerten fc MacFadveru 120111) and decollima tes. Early attempts assumed a homoge- 
neous shock front (lRhoadall999l) or spherical outfiow (IHuang et al.lll999l ). while even recently 



studies (e.g. iGranot fc Piranll201l[ ) do not account for the radial structure of the jet. 



The practical implications of the fact that the blast wave evolution is determined by a 
very small number of variables such that scalings between different explosion energies and 
circumburst density are possible were not fully realized until very recently. The scalings 
apply in the asymptotic self-similar limits, but also in the intermediate regime where the 
two-dimensional nature of jet decoUimation is in full effect. This made it possible to use only 
a small set of simulations for different initial jet opening angles as a basis for a simulati on- 
based fit code that can be applied to broadband afterglow data (jVan Eerten et al.l 1201 if ). 



But even though a complete recalculation of the dynamics of the blast wave is no longer 
necessary, there remained the equations of radiative transfer of (a representative number of) 
rays through the evolving jet, that have to be solved for each datapoint fit iteration. For a 
large number of iterations and datapoints this procedure remains computationally expensive 
and requires the use of a parallel computer. 

The current study shows that the calculation time for a given light curve or spectrum can 
be further reduced. We demonstrate that scalability between blast waves has straightforward 
implications for scalability between light curves. In ^ we describe how the scaling relations 
for the dynamics of blast waves can be used to scale between light curves as well. We 
show that the scalings remain unchanged between the BM and ST regimes, and in ^ we 
demonstrate numerically that the scalings also hold in the intermediate regime and for off- 
axis observers. We discuss our findings in ^ 
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2. Full evolution scaling relations 

We take the BM solution for the impulsive injection of an energy Eiso in a homogeneous 
medium with density po = norrip (with uq the number density and nip the proton mass) as 
the intitial condition. In this paper we will discuss only adiabatic blast waves exploding in 
a homogeneous interstellar medium, but all conclusions drawn here can be generalized to a 
stellar wind environment as well, which will be presented in a forthcoming publication. Only 
a small number of independent dimensionless combinations of the variables determining the 
fluid state at distance r, angle 6 and source frame time te exists: A = r/cte (with c the speed 
of light), B = Eisotl/por^, and 6^0 (the initial jet half-opening angle). Any fluid quantity 
can be expressed as a dimensionless combination (e.g. n/riQ for the local number density 
n) and therefore as a function of the dimensionless variables. It follows that dimensionless 
fluid quantities are inva riant under scalings that leave the dimensionless variables invariant. 



Van Eerten et al.l (1201 ll ) made practical use of the scalings 
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The synchrotron emission spectra from the expanding blast wave can be described lo- 
cally by a series of connecting power laws, see Fig. [H Below the critical frequency Va the 
medium is optically thick due to synchrotron self-absorption. Critical frequency z/^ marks 
the synchrotron break frequency. Above the cooling-break frequency v^. the accelerated elec- 
trons lose their energy too quickly to radiate fully at these frequencies. In detailed numerical 
models the evolution of the electron distribution is traced explicitly as it advects into th e 



Downes et al. 



non-radiating, slow-mov ing and dilute downstream region, see e.g 
Van Eerten et al.l (120 lOf ). Alternatively, a steady state for the radiating flu id is assunaed ari d 
the global cooling time tc is equated to the duration of t he explosion (ISari et al.lll998l ). 



Both approaches lead to qualitatively the same behavior ( Van Eerten et al. 201o[ ). 



Synchrotron radiation from shock-accelerated electrons in a shock-generated magnetic 
fleld is parametrized as follows: p denotes the power law slope of the shock-accelerated 
electron distribution, the fraction of magnetic energy relative to thermal energy, the 
fraction of downstream thermal energy density in the accelerated electrons, and the frac- 
tion of the downstream particle number density that participates in the shock-acceleration 
process. Typical values for these parameters are p ~ 2.5, ~ 0.01, e^. ~ 0.1, ~ 1- These 
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Fig. 1. — Synchrotron spectra with different orderings of th e critical frequencies. The dif- 
ferent power law segments are labeled by letters (following iGranot fc Saril |2002| ) and the 
spectral slopes are indicated in the plots. 
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parameters can be constrained from broadband afterglow data. At late times, maintain- 
ing (^7v ~ 1 sometimes leads to unphysically low values for the lower cut-off Lorentz factor 
of the accelerated parti cle distribution (which determines Um)- For a discussion, see e.g. 
Van Eerten et al.l (|2010[ ). which explores (scale-invariant) evolution of .^tv- 



The synchrotron emission coefficient ji, for a local distribution of particles follows ju oc 
^NnBf{h',h'rn,^c)/l'^{^ — Plj)"^ in the source frame, where number density n and magnetic 
field B are in the frame comoving with the fiuid, 7 denotes the fiuid Lorentz factor, /3 the 
fiuid velocity divided by c and /z the cosine of the angle between the fiuid velocity and 
the observer direction. The function /(i/, z/^, t'c) is given by the synchrotron spectrum as 
shown in Fig [1] (without self-absorption, which is included in the absorption coefficient when 
enabled) and is normalized to one where the spectrum peaks. In the optically thin case, the 
observed fiux is given by 

F=^-^ [ d^jW7(l-/3/i)- (2) 

Here di is the luminosity distance and z the redshift. Although self-absorption is shown 
in Fig [1] for the case where z/q < z/^ < t'c, we will concentrate on the case where observer 
frequency u > Ua and leave a detailed treatment of self-absorption to future work. The 
leading order dependencies of the observed fiux on the model parameters can be calculated 
from eq. [2], the shape of the synchrotron spectrum and a model for the dynamics that sets 
the radius of the afterglow blast wave as well as the post-shock val ues of n and -B, and 



Uc (the latter two via their dependence on the local fiuid state, e.g. ISari et al.lll998l ). These 
are shown in table [H 

So far, no new arguments have been presented. However, eq. IB directly implies that 
scale invariance on the dynamical level leads to scale invariance for the flux within a given 
spectral regime (e.g. when ji, depends on z/^, i^c and z/ according to a single combination of 
power laws). The table shows this explicitly for the various spectral regimes and the BM and 
ST limits. Additionally, we have expressed t and u at z = instead of the observer frame, 
which emphasizes that different redshift scalings of the fiuxes between the different spectral 
regimes are a feature of the frame in which frequency and time are expressed. Light curves 
and spectra can be mapped onto different redshifts in a straightforward manner via a single 
{1 + z) dependency. Finally, also from eq. [21 the dependence on the radiation parameters 
eg, and remains unchanged at all times. 

The fact that observer time t scales the same as te and r (that have been integrated over 
in order to obtain the flux) follows from the t = t^ — fir /c constraint that matches different 
emission times to a single arrival time. In table [1] we have tabulated the effect of applying 
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Table 1: Scalings for flux in different spectral regimes, both in the relativistic BM limit and 
non-relativistic ST limit. Note that t and v are expressed in the frame where redshift 2 = 0. 
In the observer frame, = z//(l + z) and = t{l + z). The k and A columns denote the 
corresponding scaling of the flux under a scaling of energy, time and density according to 
eq. [31 The scaling of Ua applies to the case where Ua < < ^c- 
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the mapping 



^iso ~ '^Eisoi 



Uq = Ano, 

t' = {K/Xy/H, (3) 

to the observed flux. This scahng remains unchanged between early time BM and late time 
ST for all spectral regimes, e.g. F^^bm/ ^d,bm = Pd,stI Fd,st = kA^/^, independent of jet 
opening angle or observer angle. 



3. Numerical verification 



Fig. [2] demonstrates that the scaling between energies applies at all times. The light 
curves remain in spectral regime G throughout the ir evolution, so F'(E'-^„,t') = KF( Ei,^n,t). 
All light curves were taken from the dataset used in lVan Eerten fc MacFadyenl (j201l[ ). avail- 



able from the on-line afterglow library at http://cosmo.nyu.edu/afterglowlibrary 



To generate light curves that do not lie on the asymptote of a single spectral regime 
we need to take the spectral evolution into account. As shown in table [H both the critical 
frequencies and peak fluxes obey the same scalings throughout the blast wave evolution. 
Therefore, once we know their time evolution for a given set of parameters {Eiso, tiq, 9q) 
we can use these as a baseline to generate the connected power-law spectra for arbitrary 
values of the jet parameters. In Fig. Elwe plot the time evolution of Fp^aki ^c, for 9o = 0.2 
rad and an on- ax i s obs erver. They have been calculated using the methods presented in 
Van Eerten et al.l (l201l[ ) and are ultimately based on the simulations described therein. All 
simulations were performed using the RAM adaptive-me sh refinement ( AMR) relativistisc 



hydro dynamics (RHD) code (jZhang fc MacFadyenl |2006| ). For details, see IVan Eerten et al. 

feoiih . 



At early times i>c and Fpeak differ numerically from their asymptotically expected time 
evolution because the blast wave is initially underresolved in the RHD simulations. However, 
all explosion energy is included in the initial conditions of the simul ations, and the resulting 



drop in blast wave Lorentz factor is temporary (for more details, see lVan Eerten et al.ll2011 



Zhang fc MacFadyenl 120091 ). But initially the observed flux level and critical frequencies 
are impacted due to the temporary decrease in beaming. However, the scalability of the jet 
dynamics justifies the computational cost of extremely high resolution blast wave simulations 
and these will be presented in future work. 



For jet half-opening angle 6q = 0.2, a jet break occurs around ~ 10 days. As a result. 
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Fig. 2. — Demonstrations of scaling between optical light curves (z/ = 4.56 x lO-*^^ Hz) with 
Eiso = 10^® erg and Eiso — 10^° erg, for different observer angles 9obs- Other parameters are 
set as follows: no = 10"^ cm"^, 9q = 0.2 rad, p = 2.5, = 0.01, eg = 0.1, — 1-0, z — 0, 
= 10^^ cm. The legend in the top plot refers to all plots. 
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Fig. 3. — Evolution of critical quantities for on-axis observers. The following parameters 
were used: Eiso = 5 x 10^^ erg, no = 10~^ cm"^, = 0.2 rad, p = 2.5, = 0.01, eg = 0.1, 
^jv = 1-0, z = 0, cIl = 10^^ cm. From top to bottom, the plots show z/^, Uc and Fp^ak- The 
BM and ST asymptotic slopes are indicated by dashed lines. The dotted li ne in the center of 
the Uc plot indicates the temporal slope of Uc found by iFilgas et al.l (l201l[ l for GRB 091127. 
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the time evolutions change significantly. The temporal slope of z/^ turns over quickly from 
u„i oc (BM) to the far steeper ST slope z/^ oc t~^. The temporal slope for Uc is not 

only less steep at late times, its late time ST asymptote also lies significantly higher than 
the early time BM asymptote, and as a result Uc actually rises for some time after the 

—3/5 

jet break. This effect will be less severe for larger opening angles, since Uc oc E^^J in 

—3/5 

the BM regime and Uc oc Ej in the ST regime, where Ej the total energy in both jets 
(and therefore in the final ST sphere). The two energies are related via Ej ^ Eisod'o/'^- 
It is also worth noting that, before rising, the temporal slope of temporarily steepens 
beyond —1/2. A steepening of the cooling break frequency to Uc oc has recentl y been 



observed in GRB 091127 by comparing optical and X-ray data (iFilgas et al.ll201l[ ). Our 
plot shows that this is, in principle, not inconsistent with simulations (and therefore with 
the standard model, since we do not expand upon the standard synchrotron framework by 
including features like evolving microphysics parameters such as e^). However, we caution 
against overinter pretation of the post-break z/c evolution because our approach to electron 



cooling (based on lSari et al.lll998l ) relies on a single global cooling time approximation rather 



than on tracing the local accelerated ele ctron distribution (for a comparison between the two 



approaches, see IVan Eerten et al.ll2010l . In the example there, z/^ for local cooling is typically 
higher by a factor ~ 5). Given this caveat for Uc, a clear steepening of and immediately 
post jet break is a general prediction of our study, with the steepening of Um being more 
robust. 

The final feature in all three evolution plots is the onset of the counterjet around ~ 250 
days, resulting in a relative increase of Fp^ak and z/^ and a decrease in z/^. This effect is 
strongest around ~ 1500 days. 

Using both the on-axis baselines shown in Fig[3]and the baselines for 6q = 0.2 rad, 6obs = 
26q/3, we have generated spectra (excluding synchrotron self-absorption) for a different set 
of explosion and radiation parameters and compare these to spectra calculated directly from 
simulations. The results are shown in Fig |H The off-axis angle is equal to the average 
observer angle assuming randomly oriented jets and no detection if 9obs > ^o- The scaling 
approach correctly captures the peak flux and break frequencies. The scalings-based spectra 
can be further improved upon by including smooth power law transitions between different 
spectral regimes. Fig [3] suggests that some dependency on 9obs is to be expected. 



4. Summary and Discussion 



We show that gamma-ray burst afterglow spectra and light curves above the synchrotron 
self-absorption break can be generated for arbitrary explosion and radiation parameters by 
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Fig. 4. — Comparison between simulation and scalings based spectra. The following param- 
eters are used for the spectra: Eiso = 35 x 10^^ erg, no = 11 x 10~^ cm~^, — 0.2 rad, 
p — 2.5, €b — 0.03, Ce = 0.1, ^jv = 1-0, z = 0, cIl = 10^* cm. The top plot shows a spectrum 
viewed on-axis, the bottom plot shows the spectrum observed from the average expected 
observer angle. The spectra are taken at an observer time t = 43 days. 



- 12 - 



scaling the values of a few key parameters [Fpeak, J^c, ^m) from a given baseline. The baseline 
only needs to be calculated once for each observer and jet coUimation angle. In the current 
study we have used sharp transitions between the differe nt spectral regimes, but smooth 
power law transitions can be parametrized (e.g. following iGranot &: Sarill2002l ). Although 
we have confined our study to a homogeneous circumburst environment, the generalization 
to a stellar wind environment is straightforward and will be presented in future work. The 
blast waves initially follow the ultra-relativistic Blandford-McKee self-similar solution and 
gradually spread out to the late time non-relativistic Sedov-von Neumann- Taylor stage. 

We plot the time evolution of the key parameters Fpeak, J^m, and the plots reveal that 
the critical frequencies are strongly affected by the jet break. After the jet break, the Um 
temporal slope quickly drops to the steep late time Sedov- Taylor slope z/^ oc t~^, while the 
cooling break first steepens then rises to meet the level of its late tir ne asymptote. Th e 
steepening of the temporal slope of z/c has been observed for GRB 091127 ( iFilgas et al.ll201ll ). 
though we caution against overinterpretation of our re sults given our s implified approach to 
electron cooling using a global cooling time (following ISari et al.lll998[ ). 



The scaling-based light curves fully include all two-dimensional dynamical features of 
expanding and decelerating afterglow blast waves. The fact that light curves can now be 
instantly generated for the standard synchrotron afterglow model, while taking realistic dy- 
namics into account, has the potential to strongly impact afterglow light curve fitting. In 
future work, we will implement the time evolution of key parameters for a wide range of 
observer and jet coUimation angles into a fit code that will be made available for download 
on http://cosmo.nyu.edu/afterglowlibraryj, Our findings imply that it is now possible 
in principle to fully fit for the shape of the jet break (e.g. at early time X-ray and optical 
frequencies, see lEvans et al.ll2009l : iRacusin et al.ll2009l for examples from S wift). More gener- 



ally, t he accuracy of differ ent parametrizations of the jet break transition ( iBeuermann et al. 



19991 : iHarrison et al.lll999l ) can be assessed. With the general shape of the transition into the 
ST regime being known for arbitrary explosion parameters, the exact moment when the ST 
asymptote is reached in the light curve (see e.g. iLivio fc Waxmanll2000l : IWveoda et al.ll201l|) 



is no longer relevant, which w ill improve late time radio calorimetry (see e.g. iBerger et al. 
20041 : IShivvers fc Bergeilboilh . 



In this work we have not discussed synchrotron self-absorption in detail, but we do show 
that the critical frequency Ua obeys the same scaling relations as the other key parameters. 
The effect of self-absorption will be investigated in a follow-up study. Although the scaling 
of Ua is encouraging from the perspective of broadband afterglow fitting, features such as 
the chromaticity of the jet-break across the self-absorption break (jvan Eerten et al" 2011) 
might impact the spectral slope and the sharpness of the transition into the self-absorption 
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spectral regime, rendering a parametrization slightly more complex. 



In our previous study (jVan Eerten et al.ll201ll ) we presented a method to quickly cal- 
culate light curves based on scalings on the level of the jet dynamics. This method had the 
disadvantage that a full radiative transfer was still required for each datapoint. As a result, 
broadband afterglow fitting-based on this approach still requires use of a large parallel com- 
puter in practice. The current method no longer requires radiative transfer calculations and 
is therefore vastly superior in ter ms of computational cos t. Nevertheless, in two key aspects 
the 'BOX-based' framework from IVan Eerten et al.l (120111 ) remains relevant. First, since the 
scalings there happen on the level of the dynamics, no parametrizations are necessary in 
order to describe the transitions between different spectral regimes and smoothly connected 
power laws emerge naturally. Second, the scalable BOX-based blast wave dynamics data 
provide a testing lab for studying the effect of various radiative processes, while the current 
study takes synchrotron radiation as its starting point. 
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